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Abstract. This paper is concerned with several not-quantum aspects of black holes, 
with emphasis on theoretical and mathematical issues related to numerical modeling 
of black hole space-times. Part of the material has a review character, but some new 
results or proposals are also presented. We review the experimental evidence for ex- 
istence of black holes. We propose a definition of black hole region for any theory 
governed by a symmetric hyperbolic system of equations. Our definition reproduces 
the usual one for gravity, and leads to the one associated with the Unruh metric in 
the case of Euler equations. We review the global conditions which have been used 
in the Scri-based definition of a black hole and point out the deficiencies of the Scri 
approach. Various results on the structure of horizons and apparent horizons are pre- 
sented, and a new proof of semi-convexity of horizons based on a variational principle 
is given. Recent results on the classification of stationary singularity- free vacuum so- 
lutions are reviewed. Two new frameworks for discussing black holes are proposed: a 
"naive approach", based on coordinate systems, and a "quasi- local approach", based 
on timelike boundaries satisfying a null convexity condition. Some properties of the 
resulting black holes are established, including an area theorem, topology theorems, 
and an approximation theorem for the location of the horizon. 



1.1 Introduction 

Black holes belong to the most fascinating objects predicted by Einstein's theory 
of gravitation. Although they have been studied for years,^ they still attract 
quite a lot of attention in the physics and astrophysics literature and, in fact, an 
exponential growth of the number of related papers can be observed.^ It turns 
out that several field theories are known to possess solutions which exhibit black 
hole properties: 

• The "standard" gravitational ones which, according to our current postu- 
lates, are black holes for all classical fields. 

• The "dumb holes", which are the sonic counterparts of black holes, first 



discussed by Unruh |127| 



* Partially supported by a Polish Research Council grant # KBN 2 P03B 130 16. 
^ The reader is referred to the introduction to |^ for an excellent concise review of 

the history of the concept of a black hole, and to [^,|8^ for more detailed ones. 
^ A search on the key word "black holes" on 2.XI.2001 reveals 213 papers on gr-qc for 

the current year and 773 papers from 1991, date of the beginning of the archive; the 

figures on astro-ph and hep-th are very similar. 
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The "optical" ones - the black-hole counterparts arising in the theory of 
moving dielectric media, or in non- linear electrodynamics ^5, 111|. 
The "numerical black holes" - objects constructed by numerical general rel- 
ativists. We shall argue below that this leads to the need of introducing, and 
studying, new frameworks for the notion of black holes; two such frameworks 
("naive black holes" and "quasi- local black holes") will be introduced and 
studied in Sections 1.8.1 and 1.8.2 below. 



(An even longer list of models and submodels can be found in Q.) In this 
work we shall discuss various aspects of the above. The reader is referred to 
[^^,1?^,^, 113, 133 and references therein for a review of quantum aspects of 
black holes. Let us start with a short review of the observational status of black 
holes in astrophysics. 



1.2 Experimental evidence 

While there is growing evidence that black holes do indeed exist in astrophysical 
objects, and that alternative explanations for the observations discussed below 
seem less convincing, it should be borne in mind that no undisputed evidence 
of occurrence of black holes has been presented so far. The flagship black hole 
candidate used to be Cygnus X-1, known and studied for years (c/., e.g., |p7| ). 
Its published mass has been going up and down over time: an optimistic inter- 
pretation of this phenomenon is that there has been mass accretion during the 
ascending periods, and Hawking radiation during the descending ones; a more 
realistic one is that there is still considerable uncertainty in the determination 



of this mass. Table LP lists a series of very strong black hole candidates in 
X-ray binary systems; Mc is mass of the compact object and M^, is that of its 
optical companion; some other candidates, as well as references, can be found 



in [103, 109]. The binaries have been divided into two families: the High Mass 
X-ray Binaries (HMXB), where the companion star is of (relatively) high mass, 
and the Low Mass X-ray Binaries (LMXB), where the companion is typically 
below a solar mass. The LMXB's include the "X-ray transients", so-called be- 
cause of flaring-up behaviour. This particularity allows to make detailed studies 
of their optical properties during the quiescent periods, which would be impossi- 
ble during the periods of intense X-ray activity. The stellar systems listed have 
X-ray spectra which are neither periodic (that would correspond to a rotating 
neutron star), nor recurrent (which is interpreted as thermonuclear explosions 
on a neutron star's hard surface). The final selection criterion is that of the mass 
Mc exceeding the Chandrasekhar limit Mc ~ 3 solar masses According to 

the authors of [^, the strongest black hole candidate in 1999 was V404 Cygni, 
which belongs to the LMXB class. 



Table 1.1 should be put into perspective by realizing that, by some esti- 
mates |99||, a typical galaxy - such as ours - should harbour 10^ — 10* stellar 

^ The rece nt review [lOE] lists thirteen black hole candidates. 

* See 1 103 for a discussion and references concerning the value of Mc. 
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Type 


Binary system 


Mc/Mq 




HMXR- 




11-21 


24-42 




LMC X-3 


5.6 -7.8 


20 




LMC X-1 


> 4 


4-8 


LMXB: 


V 404 Cyg 


10-15 


^ 0.6 




A 0620-00 


5-17 


0.2-0.7 




GS 1124-68 (Nova Muse) 


4.2-6.5 


0.5-0.8 




GS 2000+25 (Nova Vul 88) 


6-14 


^ 0.7 




GRO J 1655-40 


4.5 - 6.5 


^ 1.2 




H 1705-25 (Nova Oph 77) 


5-9 


^ 0.4 




J 04224-^32 


6-14 


^ 0.3 - 0.6 



black mass holes. We note an interesting proposal, put forward in ||2^, to carry 
out observations by gravitational microlensing of some 20 000 stellar mass black 
holes that are predicted |105| ] to cluster within 0.7 pc of Sgr A* (the centre of 
our galaxy). 

It is now widely accepted that quasars and active galactic nuclei are pow- 
ered by accretion onto massive black holes |101, 135 1. Further, over the last few 
years there has been increasing evidence th at massi ve dark objects may reside 
at the centres of most, if not all, galaxies [100, 12C]. In several cases the best 
explanation for the nature of those objects is that they are "heavyweight" black 
holes, with masses ranging from 10^ to 10^° solar masses. Table 1.2^ lists some 
supermassive black hole candidates; some other candidates, as well as precise 
references, can be found in p7, 102, 104, 119|. The main criterion for finding can- 



didates for such black holes is the presence of a large mass within a small region; 
this is determined by maser line spectroscopy, gas spectroscopy, or by measuring 
the motion of stars orbiting around the galactic nucleus. The reader is referred 
to ]10^ for a discussion of the maser emission lines and their analysis for the 
supermassive black hole candidate NGC 4258. An example of measurements via 
gas spectrography is given by the analysis of the Hubble Space Telescope (HST) 
observations of the radio galaxy M 87 |126] (compare |101]): A spectral analysis 
shows the presence of a disk-like structure of ionized gas in the innermost few 
arc seconds in the vicinity of t he n ucleus of M 87. The velocity of the gas mea- 
sured by spectroscopy (cf . Fig. LI ) at a distance from the nucleus of the order of 
6 X 10^^ m, shows that the gas recedes from us on one side, and approaches us on 



^ The table lists those galaxies whi ch a re listed both in [104 
some candidates from earlier lists 



and 

19| ] do not occur any more in 
of the observations listed have been published in 2000 or 2001. 




we note that 
Nineteen 
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Table 1.2. Twenty-nine supermassive black hole candidates (from [p7, 
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) 



dynamics of 


host galaxy 


Mh/Mq 


host galaxy 


Mh/Mq 


water maser discs: 


NGC 4258 


4 X 10^ 










9 V 1 n** 


M 87 
iVi * 


^ V 1 

A lU 




TSJpp 9707 


/I V 1 n"^ 




9 V 1 

Z A lU 






V 1 

A lU 




4 V 1 

^ A lU 




isjpp CI 90 

IN VjrVy OIZO 


9 V 1 

Z A lU 


NPr* fl9Cil 

IN UZOi 


X lU 






V 1 

A lU 






stars: 


NGC 821 


4 X 10^ 


NGC 1023 


4 X 10^ 




NGC 2778 


1 X 10^ 


NGC 3115 


1 X lO'' 




NGC 3377 


1 X 10** 


NGC 3379 


1 X 10* 




NGC 3384 


1 X 10^ 


NGC 3608 


1 X 10* 




NGC 4291 


2 X 10** 


NGC 4342 


3 X 10* 




NGC 4473 


1 X 10** 


NGC 4486B 


5 X 10* 




NGC 4564 


6 X 10^ 


NGC 4649 


2 X lO'' 




NGC 4697 


2 X 10** 


NGC 4742 


1 X lO'' 




NGC 5845 


3 X 10** 


NGC 7457 


4 X 10^ 




Milky Way 


2.5 X 10'^ 







the other, with a velocity difTerence of about 920 km s ^ . This leads to a mass 
of the central object of ~ 3 x 10^ Af©, and no form of matter can occupy such 



a small region except for a black hole. Figure 1.2 shows another image, recon- 
structed out of HST observations, of a recent candidate for a supermassive black 
hole - the (active) galactic nucleus of NGC 4438 [Q. To close the discussion of 
Table |1.2| , we note that the determination of mass of the galactic nuclei via direct 
measurements of star motions has been made possible by the unprecedentedly 
high angular resolution and sensitivity of the HST. There seems to be consensus 



[p7, 104, 12c] that the two most convincing supermassive black hole candidates 



are the galactic nuclei of NGC 4258 and of our own Milky Way. 

There is a huge gap between the masses of stellar mass black holes and 
those of the supermassive black holes and, in Bromley's terminology |^0| , "mid- 
dleweight" candidates would be welcome. Such a tentative black hole with M of 



the order of 460Mq has been identified in M 82 by A. Ptak and R. Griffiths |117|] . 
E. Colbert and R. Mushotzky give a list of compact X-ray sources which 
could include some further black holes with masses in the 10^ — W^Mq range. 
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Let us close this section by pointing out the review paper p5[ | which discusses 
both theoretical and experimental issues concerning primordial black holes. 



1.3 Causality for symmetric hyperbolic systems 

The usual, Scri based, definition of a black hole proceeds as follows: one starts 
by introducing a preferred asymptotic region - this is usually taken to be a 
conformal boundary. Then the black hole region is defined as the set of points 
which cannot send information to the asymptotic region; we will return to the 



details of this construction in Section 1.4 below. The key ingredient here is 



the notion of "not being able to send information" - this is usually defined 
using the metric, together with the postulate that the propagation speed for 
objects carrying information is bounded from above by that of null geodesies. 
The aim of this section is to provide a precise mathematical version of this notion, 
via Proposition Q below, for any symmetric hyperbolic system of differential 
equations, without introducing any supplementary postulates. The point is that 
every physical theory defines its own ca usalit y via the corresponding system of 



equations. We will then show in Section 1.3.1 how the "Unruh metric" arises in 
the application of our construction to the Euler equations. 

Let, thus, L be a quasi-linear, symmetric hyperbolic, first order partial dif- 
ferential operator acting on sections (p oi a, bundle E, with scalar product (•, •), 




Fig. 1.1. Hubble Space Telescope observations of spectra of gas in t he v icinity of the 
nucleus of the radio galaxy M 87, NASA and H. Ford (STScI/JHU) ||l23|. 
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over a manifold ^\ in a local trivialization over a coordinate patch ^ ^ we 
have 

m:^A^^{ip,xP)d^^ + B{ip), (1.1) 

where the A'^'s are endomorphisms of the fibers of symmetric with respect to 
(•,•). Given a section Lp of E (defined perhaps on a subset of ./#), a hypersurface 
is said to be spaceZifc^ if one of its fields of co-normals satisfies 

this is understood in the sense of strict positive definiteness of endomorphisms. 
By definition of a symmetric hyperbolic system, every point of the domain of 
definition of lies on some spacelike hypersurface. For p € ^ we set 

Sr^ {a e | [lp (p) ,p) is positive definite} C T*J^ , (1.2a) 

"^p := =3^'' C {a e I {lp {p) ,p) is non-negative definite} C T*.^ , 

(1.2b) 

and 

:=Upe^^;cT*^, Upe^^p ^ ■ (1-3) 

The bundle is thus the bundle of covectors normal to spacelike hypersurfaces. 
Because positive definite matrices form a convex set, each set is a convex 

® In quasi-linear theories, in which the coefficients of the highest derivatives depend 
upon the fields, the notion of spacelikeness will of course depend upon the given set 
of fields. 




ve Galaxy NGC 4438 

Space Telescope ■ WFPC2 



Fig. 1.2. Hubble Space Telescope observations [g5| of the nucleus of the galaxy NGC 
4438, from the STScI Public Archive ]l23| . 
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cone in T^Ji. The bundles ^« C TJi, respectively "^f C TJi, of timehke, 
respectively causal, vectors are defined by duality: 

^ptt {X e |V a e =^p^ X^a^ > 0} , (1.4a) 

'rfi ■= ^{X e Tp^ |V a e ,^p^ X^a^ > 0} . (1.4b) 

Recall that all standard classical field equations, i.e., Euler's equations, the scalar 
wave equation, the wave map equation. Maxwell's equations, Yang-Mills equa- 
tions, or Einstein's equations, can be written as a symmetric hyperbolic system 



(cf., e.g., |TTq] or [125, Vol. IV]). For example, the wave equation for a scalar 



field It on a Lorentzian manifold {^,g) can be written in the form (1.1) by 
introducing a new set of variables 

LP := (u,e^(u)) , 

where the e^'s, /i = 0, . . . , n, form an ON- frame for g. In this case i^" coincides 
with the bundle of timelike future pointing vectors of the metric g (where the 
future direction is determined by that of eo), coincides with the bundle of 
causal future pointing or vanishing vectors of g, while and are, modulo 
sign (we are using the signature (— , . . . , -I-)), "i^" and '^'^ with indices lowered 
using the metric" . The same remains true for wave maps, for linear electrody- 
namics (considered as a first order system for E and B), for the Yang-Mills 
equations in the Lorentz gauge, and for the usual symmetric hyperbolic reduc- 
tion of Einstein's equations based on harmonic coordinates. (See the end of this 
section for more comments about the Einstein case.) 

We emphasize that the bundles JJ^", etc., are defined solely by the system of 
equations under consideration, independently of any metric. They can be used 
to determine the causality properties of a given symmetric hyperbolic system 
by mimicking the usual metric constructions ]T3|,[7S , 112. 115]: timelike future di- 



rected paths are defined as piecewise differentiable maps 7 from / to where 
/ C M is an interval, with the property that the tangent 7 is in wherever 
defined; causal future directed paths are defined as Lipschitz continuous paths 
such that 7 is in '■^'^ wherever defined. Past directed paths are paths which are ob- 
tained from future directed paths by a reversal of parameterization. The timelike 
future I'^{f2), respectively the causal future J^(f2), of a set Q C ^ is defined 
as the set of points which can be reached form fi by following a future directed 
timelike, respectively causal, path. Causal and timelike futures and /~ are 
obtained form the above definition by replacing "future" by "past" . Notions such 
as global hyperbolicity, and so on, can be defined in the usual way. It might be 
convenient to write I^{f2,L), etc., to emphasize the dependence upon the sys- 
tem of equations L, if ambiguities can arise. We note that for semi-linear systems 



- by definition, these are the systems in which the coefficients A'^ in (1.1) do not 
depend upon ip ~ the resulting causal constructs are independent of the solution 
if under consideration; however, this will not be the case in general. 

We are ready now to address the issue of main interest in this work - the 
notion of a black hole. Given a set f2 C ^ and a symmetric hyperbolic system 
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L we define the black hole region 3§n (L) of Q as 

:3§n{L) := ^\J-{n,L) . (1.5) 

The object so defined depends upon the set i7, and acquires its full meaning 
when fi is naturally distinguished by the problem at hand. 

In several cases, some of which have already been mentioned, the causality 



theory constructed above arises out of a metric; we will show in Section |1.3.1 
below that this holds also for the Euler equations. The intuitive meaning asso- 
ciated to the notion of a black hole region is that no information from 
can reach fl. This is made precise by the following result, the proof of which 
proceeds by causality arguments which are known in principle: 

Proposition 1 Let if he a solution of 

LM = J , (1-6) 



where L is symmetric hyperbolic operator (1.1) on a manifold , and suppose 
that there exists a Lorentzian metric g such that the sets defined in (1.4b) are 
future light cones of g. Assume that {^,g) is globally hyperbolic with Cauchy 
surface ^ and consider a set f2 C ^ such that il§^2{L) 7^ 0. If ip' is a solution 
of 

L[^'] = J' with (^1^ , (1-7) 

and if the difference of the source terms J ~ J' is supported in .^^2(^)1 then 
if — ip' on J" ( J7, L) . 

We believe that the conclusion above is true without the assumption that 
causality is determined by a metric, but we are not aware of any studies of this 
question in the current context; it would be of interest to settle this. The reader 
is referred to p7[ | for an analogous analysis of strictly hyperbolic systems, and 
to for that of second order Lagrangean hyperbolic systems. 

We close this section by mentioning that care has to be taken when inter- 
preting the notions above for systems of equations which are not directly in first 
order symmetric hyperbolic form: it can happen that the physical equations can 
be then rewritten in several different symmetric hyperbolic forms, leading to dif- 
ferent notions of pasts, futures, black hole regions, etc. This does indeed happen 
for the Einstein equations: in Section 4.1] a symmetric hyperbolic system is 
presented which can be used to solve the Einstein equations and in which the 
causal cones, as defined above, differ from the metric ones; see the discussion 
after Equation (4.25) there. Another example of this kind arises if the T- vector 
of the system of equations given in Sect. 5.2 of is chosen to be spacelike. 
This is related to the acausal (in the usual metric sense) propagation of gauge 
degrees of freedom in the associated systems, and does of course not affect the 
propagation of any physically relevant quantities. Given a system of equations, 
say (S), a simple way out of this problem is to define the physical black hole 
region as the intersection, as L runs over the set of those symmetric hyperbolic 
systems L the solutions of which include solutions of (S), of the associated black 
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hole regions .'^n{L). We note that the existence of exact solutions such as the 
pp- waves, where nonlinear perturbations propagate along null geodesies, guaran- 
tees that every symmetric hyperbolic system L which reproduces the Einstein's 
equations will necessary contain the metric light cone as its causal light cone 
associated to L. It follows that the above prescription does indeed reproduce the 
usual metric definition of a black hole region for the vacuum Einstein equations. 



1.3.1 Dumb holes 



Let us illustrate the considerations of the previous section in the case of the 
Euler equations 



p (v + V • Vv) + Vp[p) = 
/) + V(pv) = . 



(1.8) 



For strictly posi tive p's (1^) can be rewritten as a symmetric hyperbolic system 
of the form (lA) with B = 0, for a field (p — (f % p), by introducing the matrices 



pidR3 



p'5] ' 



(1.9) 



which are clearly symmetric with respect to the standard scalar product on 
X R 9 {v,p)- Straightforward algebra shows that A'^'-n^ is strictly positive if 
and only if p > 0, dp/ dp > and 



{nn+v'n,f +p'^{n,f < 



g'^'Uf^n^ < , 



(1.10) 



where the Unruh metric g is defined as |127, 13 



CsP 



-I -v' 



p 

9111^ = — 

Cs 



dp 



(1.11) 



" dp ■ 

This metric exhibits a typical black hole behaviour for stationary solutions in 
which the speed v of the fluid velocity meets the speed of sound across a 
smooth hypersurface. 

The original argument of Unruh leading to (1.11) was a perturbational one, 
for irrotational flows: 

v^V^. (1.12) 



Small perturbations (j) = S<P of (p satisfy then a scalar wave equation Il2 



-gg^^d^c^) = 



in the metric (1.11). Our discussion above shows that neither the irrotationality 
of the flow, nor the perturbative character of the argument are essential for the 
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problem at hand. The fact that eausahty for the fuU non-hnear Euler's equations 
is governed by a metric can essentiahy be found in Courant and Hilbert, and is 
certainly well known to some researchers (H. Friedrich, private communication). 

There was renewed interest in the above because of the recent experiments 
with the Bose-Einstein condensates. Recall that in the Madelung formulation, 
the Gross-Pitaevskii equation governing the dynamics of the condensates can be 
rewritten in Euler form. There have been suggestions that some of the effects 
associated to black-hole type causality, including the analogue of the Hawking 
radiation, could be observed in such systems, see [|l^,|82 , 131| and references 



therein. While the resulting black holes are still referred to as sonic, this is rather 
misleading as the underlying Gross-Pitaevskii equation does not have anything 
to do with sound propagation in liquid or gaseous physical media. 

1.3.2 Optical holes 



According to M. Visser |12S], in a dielectric fluid "the refractive index n, the 
fluid velocity v, and the background Minkowski metric rj can be combined alge- 
braically to provide an effective metric 

This leads to black-hole effective geometries for refractive indices which exceed 1 
in some regions. Phenomenological models for this behaviour have been proposed 
by U. Leonhardt and P. Piwnicki, cf. |^ and references therein]]. Similarly, a 
(different) effective metric is obtained in nonlinear electrodynamics, leading - 



according to - to models which sometimes contain closed timelike curves. 
The status of these "analogous models" seems to be somewhat less clear than 
that of the sonic ones. 



1.3.3 Trapped surfaces 



We shall close this section by noting the discussion in j9£, 129| concerning the 



notion of trapped surfaces]^ for black hole geometries - such surfaces are useful 
detectors of black hole regions in GR with matter fields carying positive energy. 
It should be borne in mind that while causality concepts can be introduced for 
any symmetric hyperbolic system, the existence of a trapped surface signals the 
presence of a black hole region only when causality is governed by a metric satis- 
fying an energy positivity condition (and, i f in a Scri context, when Scri satisfies 



various regularity conditions, see Section lA below). In the non-gravitational 
black hole models the metric is the one arising out of the causality structure of 
the theory, and there are no reasons in general to believe that its Ricci tensor 
should have any preferred properties; to start with, it is defined only up to a 
conformal factor anyway. In addition, the " trapped surface" terminology is used 



See also URL bittp: //www . st-and. ac .uk/~www_pa/group/quantumopticE| 
® Cf. Section |l.6| for the difFerential-geometric definition of this notion 
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in ]129| for objects which are completely unrelated to the usual differential ge- 
ometric context, which is extremely misleading and obscures the problems at 
hand. 



1.4 Standard black holes 

The standard way of defining black holes is by using conformal completions: A 
pair g) is called a conformal completion of g) if ^ is a manifold with 
boundary such that: 

1. ^ is the interior of 

2. there exists a function with the property that the metric g, defined to be 
il^g on extends by continuity to the boundary of with the extended 
metric still being non-degenerate throughout, 

3. f2 is positive on differentiable on vanishes on with df2 nowhere 
vanishing on y. 

We emphasize that no assumptions about the causal nature of Scri are made 
so far. The boundary of ^ will be called Scri, denoted J^. In the usual text- 



books 1 73, 132 1 smoothness of both the conformal completion and the metric g 



is imposed, though this can be weakened for many purposes. We set 
Assumi ng var ious global regularity conditions on ^ , to which we shall return 



in Sect. 1.4.1 , one then defines the black hole region ^ as 

^:=^\J~(^+). (1.13) 

Let us point out some drawbacks of this approach: 

• Non-equivalent Scri's: Conformal completions at null infinity do not have 
to be unique, an example can be constructed as follows: the Taub-NUT 



metrics can be locally written in the form |106| 



-U^^de + {2LfUal + [e + L^){a^ + a^) , (1.14) 
[/(t) = _l+ (1.15) 

where (Ti, ct2 and are left invariant one-forms on SU{2) w S^. The con- 
stants L and m are real numbers with L > 0. Parameterizing with Euler 
angles {fi,9,ip) one is led to the following form of the metric 

g = ^U-^de + {2LfU{d^i + cosOdiff + {t^ + L^){d9^ + sin^ Odip^) . 

The standard way of performing extensions across the Cauchy horizons t± := 
M ± V Af 2 + ^2 jg introduce new coordinates 

{t, ^l, e, ^)^{t,fi± f [2LU{s)]-'ds, 9, ^) , (1.16) 

Jto 
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which gives 

g± — ±4L{dfi + cosddLp)dt 

+ {2LfU{dfi + cosedipf + [t^ + L^)(d6|2 + sin^ OdLp^) . (1.17) 

Each of the metrics g± can be smoothly conformally extended to the bound- 
ary at infinity = oo" by introducing 

x = l/t , 

so that ( l.lTp becomes 



g± = X ^(^T 4:L{dfi + cos 9dip)dx 

+ {2L)^x^U{dfi + cos9dip)^ + (1 + L^x^){d9^ + sin^ 9dip^)j (1.18) 

In each case this leads to a Scri difleomorphic to S^. There is a simple 
isometry between g_|- and g_ given by 

{x,fj,,9,(f) {x,-fi,9,-ip) 

(this does correspond to a smooth map of the region t £ (t+, oo) into itself, 
cf. 1^^), so that the two Scri's so obtained are isometric. However, in addition 
to the two ways of attaching Scri to the region t e {t^,oo) there are the 
two corresponding ways of extending this region across the Cauchy horizon 
t — t^, leading to four possible manifolds with boundary. It can then be seen, 
using e.g. the arguments of that the four possible manifolds split into 
two pairs, each of the manifolds from one pair not being isometric to one from 
the other. Taking into account the corresponding completion at "t = — oo", 
and the two extensions across the Cauchy horizon t — t^, one is led to 
four inequivalent conformal completions of each of the two inequivalent p2| 
time-oriented, maximally extended, standard Taub-NUT space-times. 

Under mild completeness conditions in the spirit of ||66|] , uniqueness of 
^+ as a point set in past- distinguishing space-times should follow from the 
TIP and TIF construction of however, uniqueness of Scri's differentiable 
structure within that framework is far from being clear. 

Yet another approach to uniqueness has been proposed by Geroch p5[ |: 
A completion as defined at the beginning of this section is said to be a 
regular asymptote if, given any point p € and any non-zero null vector 
£ e Tq^^ such that the maximally extended null geodesic with tangent £ 
does not meet p, then there are neighbourhoods ^ C ^ of p and 'V C 
of £ such that no maximally extended null geodesic with tangent in Y meets 
. This regularity condition is not satisfied by the Taub-NUT completions 
described above. Supposing that the set of regular asymptotes is non-empty, 
Geroch [jssl . Theorem 2, p. 14] gives an argument for existence of a maximal 
regular asymptote, unique up to well behaved conformal transformations. 
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Now, we have not been able to fill in the details]^ of sonic of the arguments 
suggested in ||6^, but the construction given is plausible enough so that we 
expect that it might actually be correct. Whatever the case, Geroch's regu- 
larity requirement is a global condition which seems to be difficult to control 
in general situations.]^ In particular it could happen that many space-times 
of interest admitting conformal completions do not admit any regular ones. 
It should be borne in mind that one of the main applications of the ^ frame- 
work is the possibility of uniquely defining global charges such as mass, angu- 
lar momentum, etc. (cf. |4^ ] for a new approach to this question), in which 
issues about global behaviour of all geodesies seem completely irrelevant. 
The fact that there is no suitable uniqueness result for the differentiable 
structure of conformal completions is a serious gap in our understanding of 
objects such as the Trautman-Bondi mass. 

We note that uniqueness of a class of Riemannian conformal completions 
at infinity has been established in Section 6] ; this result can probably be 
used to obtain uniqueness of differentiable structure of Lorentzian conformal 
completions for Scri's admitting cross-sections. Further partial results on the 



problem at hand can be found in [121| 



Poorly differentiable Scri's: In all standard treatments [^,|73, 132 1 it is 



assumed that both the conformal completion ^ = ^V}J' and the extended 



metric g are smooth, or have a high degree of differentiability |114|. This 
is a restriction which excludes most space-times which are asymptotically 
Minkowskian in lightlike directions, see |^,|9^] and references therein. Poor 
differentiability properties of J" change the peeling properties of the gravi- 
tational field , but most - if not all - essential properties of black holes 
should be unaffected by conformal completions with, e.g., polyhomogeneous 
differentiability properties as considered in [^Q. It should, however, be 
borne in mind that the hypothesis of smoothness has been done in the stan- 
dard treatments, so that in a complete theory the validity of various claims 
should be reexamined. Some new results concerning existence of space-times 
with a poorly differentiable Scri can be found in 

The structure of «+: The current theory of black holes is entirely based 
on intuitions originating in the Kerr and Schwarzschild geometries. In those 
space-times we have a family of preferred "stationary" observers which follow 
the orbits of the Killing vector field dt in the asymptotic region, and their 
past coincides with that of J^"*". It is customary to denote by i+ the set con- 



It is not completely clear that the map defined by Geroch in his proof is simultane- 
ously differentiable for all regular asymptotes, as asserted in p. 14]. The regularity 
condition guarantees Hausdorfi'ness of the maximal regular asymptote constructed 
in Jbs] , p. 14], but does not seem to guarantee its differentiability in any obvious way. 
Globally hyperbolic conformal completions (in the sense of manifolds with bound- 
ary) should satisfy Geroch's condition. Nevertheless, it should be borne in mind that 
good causal properties of a space-time might fail to survive the process of adding 
a conformal boundary. For example, adding a Cauchy horizon to a maximal glob- 
ally hyperbolic space-time sometimes leads to space-times which are not globally 
hyperbolic in the sense of manifolds with boundary. 



14 P. T. Chrusciel 



sisting of the points t ^ oo, where t is the KiUing time parameter for those 



observers. Now, the usual conformal diagrams for those space-times [73, 107| 
leave the highly misleading impression that i"*" is a regular point in the con- 
formally rescaled manifold, which, to the best of our knowledge, is not the 
case. In dynamical cases the situation is likely to become worse. For exam- 
ple, one can imagine black hole space-times with a conformal diagram which, 
to the future of a Cauchy hypersurface t — 0, looks as in Fig. 1.3. In that 



singularity 




Fig. 1.3. An asymptotically flat space-time with an unusual i"* 



diagram the set i+ should be thought of as the addition to the space-time 
manifold ^ of a set of points "{t = oo,q G ^}" , where t e [0, oo) is the 
proper time for a family of observers ^. The part of the boundary of ^ 
corresponding to is a singularity of the conformally rescaled metric, but 
we assume that it does not correspond to singular behaviour in the physical 
space-time. In this space-time there is the usual event horizon Si correspond- 
ing to the boundary of the past of , which is completely irrelevant for 
the family of observers ff, and an event horizon #2 which is the boundary of 
the true black hole region for the family i?, i.e., the region that is not acces- 
sible to observations for the family ff. Clearly the usual black hole definition 
carries no physically interesting information in such a setting. 
Causal regularity of Scri: As already pointed out, in order to be able 
to prove interesting results the definition ( 1.13| ) should be complemented by 
causal c onditi ons on The various approaches to this question, discussed 
in Sect. 1.4.1 , are aesthetically highly unsatisfactory: it appears reasonable 
to impose causal regularity conditions on a space-time, but why should some 
unphysical completion have any such properties? Clearly, the physical prop- 
erties of a black hole should not depend upon the causal regularity - or 
lack thereof - of some artificial boundary which is being attached to the 
space-time. While it seems reasonable and justified to restrict attention to 
space-times which possess good causal properties, it is not clear why the 
addition of artificial boundaries should preserve those properties, or even be 
consistent with them. Physically motivated restrictions are relevant when 
dealing with physical objects, they are not when non-physical constructs are 
considered. 
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• Inadequacy for numerical purposes: Several]^ numerical studies of black 
holes have been performed on numerical grids which cover finite space-time 
19, 86|. Clearly, it would be convenient to have a set-up which 



regions 

is more compatible in spirit with such calculations than the Scri one. 



We will suggest, in Sections 1.8.1 and 1.8.2 below, two approaches in which 



the above listed problems are avoided. Before doing this, let us complete the 
presentation of the usual, Scri based, approach to black holes. 



1.4.1 Scri regularity conditions, and the area theorem 



It is easily seen that the definition ( 1.13 ) is not very useful without some com- 



pleteness conditions on . For instance, denote by ^ the standard conformal 
completion!^ of the Minkowski space-time M^^'^, with ^+ w R x S^, where the 
K factor is parameterized by a Bondi coordinate u, with u — t — r on M.^'^. One 
can obtain a new completion ./#2 by replacing by an open subset thereof, 
e.g. restricting the range of u's to (— oo, 0). This will give a non-empty black hole 
region ^ in M^'^, 

^= J+(0) , 

which is clearly a completely uninteresting statement. A way out of this problem 
(as well as of some other related ones, discussed in 66 1) has been proposed!^ 



in | |66| ] , for space-times which are vacuum near the conformal boundary assuming 
that the cosmological constant is zero: R. Geroch and G. Horowitz introduce a 
preferred family of conformal factors SI such that the matrix Hessi7 of second 
covariant derivatives of i? vanishes at It is then required that both ^+ and 
y~ be diffeomorphic to K x 5'^, with the R factor corresponding to integral 
curves of which are assumed to be complete. 

Let us pass to a presentation of the causal regularity conditions which are 
imposed in Hawking and Ellis's approach [Q. Now, a set of conditions is only 
useful insofar as it allows to prove something. Since the "area theorem"|^ is 
one of the landmark theorems in the theory of black holes, we shall concentrate 
on the set of conditions from Hawking and Ellis as used in their treatment of 
the area theorem; we follow the presentation of [pS] , Appendix B]. One of the 
conditions of the Hawking-EUis area theorem |73|, Proposition 9.2.7, p. 318] is 
that spacetime {^,g) is weakly asymptotically simple and empty ("WASE", 
p. 225]). This means that there exists an open set C ^ which is isometric 

The numerical simulations in [p^, [l^,[7o| cover regions extending all the way to in- 
finity, wit hin fra meworks which seem to be closely related to the "naive" framework 
of Section 1.8.1 below. 

Inroduce coordinates u — t — r, x — 1/r, 6, cp on \ {r = 0}, with (r, 6, ip) - the 
usual spherical coordinates on R^, and set ~ {x — 0} in those coordinates. 
" The "WASE setup" of presented below, does not solve the problem, cf. |6(1|. 
Recall that the area theorem is the statement that the area of cross-sections of event 
horizons is non-increasing towards the future. One possible precise version of this 
theorem can be found in Theorem H below. 



16 P. T. Chrusciel 



to D , where is a. neighborhood of null infinity in an asymptotically 
simple and empty (ASE) spacetime {^',g') | |73| , p. 222]. ft is fm-ther assumed 
that ^ admits a partial Cauchy surface with respect to which ^ is future 
asymptotically predictable ([Q, p. 310). This is defined by the requirement that 
J^'^ is contained in the closure of the future domain of dependence I?"*" (=5^;^) 
of , where the closure is taken in the conformally completed manifold ^ = 
^ U U 1 with both and J'^ being null hypersurfaces. Next, one 
says that {^,g) is strongly future asymptotically predictable ( |f73|| , p. 313) if it 
is future asymptotically predictable and if J+(^) n J^(^+;^) is contained 
in Finally (j?^, p. 318), {^,g) is said to be a regular predictable 

space if g) is strongly future asymptotically predictable and if the following 
three conditions hold: 

(a) ,5^{^J^ {■^^', is homeomorphic to R'^\(an open set with compact closure). 
{(3) ,5^ is simply connected. 

(7) the family of hypersurfaces S^{t) constructed in |7^, Proposition 9.2.3, p. 313] 
has the property that for sufficiently large r the sets ^(r) n J~(j^"'";^) 
are contained in J+(j^~;^). 

It should be clear at this stage that this set of conditions is rather intricate. 
As it turns out, it is not clear whether or not it does suffice for a proof of 
the area theorem, as asserted in j?^. Proposition 9.2.7, p. 318]: indeed, in the 
proof of |73, Proposition 9.2.1, p. 311] (which is one of the results used in the 
proof of 1 73, Proposition 9.2.7, p. 318]) Hawking and Ellis write: "This shows 
that if W is any compact set of every generator of leaves J+(#';^)." 
The justification of this given in |7^ ] is wrong. If one is willing to impose the 
sentence in quotation marks as yet one more regularity hypothesis on then 
the arguments given in fz^ apply to prove the area theorem for black holes with 
a piecewise smooth event horizon; the results in p8|] show that this remains true 
with no supplementary conditions on the differentiability of the horizon. 

Not only is the above set of hypotheses aesthetically unappealing, it is far 
from being unique. An alternative way to guarantee that the area theorem will 



hold in the "future asymptotically predictable WASE" set-up of |73| is to impose 
some mild additional conditions on and . There are quite a few possibilities, 
one such set of conditions has been described in Appendix B]: Let ip : ^ 
'W n denote the isometry arising in the definition of the WASE spacetime 
^ . First, one requires that ■0 can be extended by continuity to a continuous 
map, still denoted by defined on . Next, one demands that there exists a 
compact set K C such that, 

V'(J+(^;^)n(9^) c J+(i^;^') , (1-19) 



see Fig. 1.4. Under those conditions the area theorem does again hold. 



Several other proposals how to modify the WASE conditions of [|73[ to obtain 



sufficient control of the space-times at hand have been proposed in ||4^, g9, IIC]. 

A completely different, and considerably simpler, treatment of the question of 
regularity of Scri needed for the area theorem, has been proposed in [p8| : A point 
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g in a set A C -B is said to be a past point of A with respect to B ii J~ {q; B)nA = 
{q}. One then says that is J^~regular if there exists a neighborhood O of 
Jif such that for every compact set C C O satisfying I~^{C; M) n ^+ 7^ there 
exists a past point with respect to J'^ in dI^{C; M) n J^"*". The condition of 
J^-regularity of is precisely what is needed for the arguments of the area 
theorem to go through. This condition is somewhat related to the i'^ -avoidance 
condition used by Galloway and Woolgar [g^ ]. 

Yet another approach has been advocated by Wald |132|: there one considers 
glo bally hyperbolic completions in a set-up that includes r\ While the conditions 
of 1 132] do lead to a coherent set-up for the validity of the area theorem, they 
introduce several new difhculties related to the low differentiability of the con- 
formally rescaled metric at i", not to mention the question of mere existence of 
such completions. This can actually be relaxed to the related weaker requirement 
of global hyperbolicity of ^ U J^'^ (in the sense of a manifold with boundary) 
which enforces J^-regularity of 0/+ |38], so that the area theorem again holds. 



1.5 Horizons 

A key notion related to the concept of a black hole is that of an event horizon, 

S -=8^ . (1.20) 

(Actually, when is not connected some care is required for a useful defini- 
tion, but this does not affect the discussion that follows.) Event horizons are 
a special case of a family of objects called future horizons: by definition, these 
are closed topological hypersurfaces threaded by null geodesies, called gen- 
erators, with no future end points, and possibly with past end points. At the 
latter, differentiability of breaks down in general; a necessary and sufficient 
condition for this breakdown has been given in jlj]. It seems that most authors 
have been taking for granted that horizons are nice, piecewise smooth hypersur- 
faces. This is, however, not the case, and examples of nowhere horizons have 
been constructed in Further, nowhere horizons are generic in the class of 
convex horizons in Minkowski space-time [ p3| . This leads to various difficulties 
when trying to study their structure, e.g., attempting to prove results such as 
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the area theorem discussed m Sect. 1.4.1, compare Theorem || below. Horizons 
are always locally Lipschitz- continuous topological hypersurfaces ]115| (compare 
Corollary [A.3| below), i.e., they are locally graphs of functions satisfying 

\Hx)-(biy)\<C\x^y\ . 

In spite of the potential low differentiability, the usual optical scalars can be 
defined almost everywhere in an Alexandrov sense, as follows: A function is said 
to be semi-convex if it is the sum of a function and of a convex function. 
It has been shown in that future horizons are semi-convex. We give an 

alternative proof of this result in an appendix: the new argument has some 
interest on its own, as it proceeds via a variational principle for horizons, which 
could be useful for numerical applications. The interest of the semi- convexity 
property relies on the following theorem of Alexandrov: 

Theorem 2 (Alexandrov |52, Appendix E]) Semi-convex functions f : 
B are "twice- differentiable" almost everywhere in the following sense: 

there exists a set Bm with full measure in B such that 

y X e Bm ^ Q e (RP)* (E) (Wy such that y ye B 

f{y) - f{x) - d.f{x){y - x) = Q{x - y,x - y) + r2{x,y) , 

with r2{x, y) = o{\x—y\'^). The symmetric quadratic form Q above will he denoted 
by ^D'^f{x), and will be called the second Alexandrov derivative of f at x. 

Points in Bm will be referred to as Alexandrov points of B. 

One can now use Theorem ^ to define the equivalent to the usual divergence 
9 of the horizon, by writing locally as the graph of a function /, and using 
the second Alexandrov derivatives of / to define (almost everywhere on 
the divergence 9m of the horizon. More precisely, let p = {t = f{q),q) be an 
Alexandrov point of J^f and let e^, i = 1, 2, 3 a basis of TpJ^ such that 

63 = K^{p)dx^^ = -dt + df{q) , g(ei,ei) = (e2, 62) = 1 , 62) = . 

One then sets 

V.K, = D^f - r^K^ , 9m = {e\e{ + e\ei)V,K, . 

where the Hessian matrix D^^f above is that of second Alexandrov derivatives. 
This reduces to the standard definition of 9 when f e C^. 

Several results of causality theory go through with 9 replaced by 9m , though 
the standard arguments sometimes have to be replaced by completely different 
ones, the reader is referred to |Q for details. In particular the Raychaudhuri 
equation, as well as the remaining optical equations, hold on almost all generators 
of Jf". 

Somewhat surprisingly, event horizons in smooth, stationary, globally hyper- 
bolic, asymptotically flat space-times satisfying the null energy condition are 
always smooth null hypersurfaces, analytic if the metric is analytic; this is a 
corollary of the area theorem of . 

Some further results concerning differentiability properties of horizons can be found 
in 
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1.6 Apparent horizons 

In spite of their name, apparent horizons are not horizons. They are usuaUy ff^ 
defined on spacehke hypersurfaces C as foUows: let H C be the set 
covered by (future) trapped surfaces: by definition, those are compact, bound- 
ary less, smooth surfaces S C with the property that 

6(3) ~X- (g'^ - n'n^)K,j < , (1.21) 

where A is the (outwards) mean extrinsic curvature of S in J/', is the outwards 
pointing unit normal to S in .y, and Kij is the extrinsic curvature of in 
There is no reason for Q to be nonempty in general; on the other hand, in 
appropriately censored space-times, a non-empty Q implies the existence of a 
black hole region. Hawking and Ellis define the apparent horizon £/ as 

£/:^df2, (1.22) 

and argue that 



= 0. (1.23) 

Their argument is correct if one assumes that jz/ is C^. However, £/ could be 
a priori a very rough set, with 9 not even defined, in which case the arguments 



of Q do not apply. Probably the simplest way out is to use (1.23) as the defi- 
nition of apparent horizon, and forge t about i7. The existence of a compact. 



boundaryless surface satisfying ( 1.2S ) does again imply the existence of a black 
hole in a Scri framework, assuming appr opriate causal properties of 

We note that if Kij = 0, then ( 1.22 ) is the equation for a minimal surface. 



Now, in the course of their proof of the Penrose inequality, G. Huisken and 
T. Ilmanen prove that the outermost minimal (in the sense of calculus of 



variations) surface is always smooth, which supports the validity of (1.23) at 
least for the outermost component of £/. However, they provide examples which 
show that this outermost minimal surface might not be the boundary of f2 in 



general. This illustrates another deficiency of the definition ( 1.22 ). Some partial 
results concerning the differentiability properties of have been obtained by 
M. Kriele and S. Hayward in KS]. R. Howard and J. Fu have recently shown 



that df2 satisfies ( 1.22| ) in a viscosity sense - this is defined as follows: Let D be 



an open set, and let p G dD. Then U is an inner (outer) support domain if U 
is an open subset of D (of CD), dU is a hypersurface and p S dU . One says 
that 9{dD) < (respectively 9{dD) > 0) in the viscosity sense if for any point 
p G dD and for any inner support (outer support) domain we have 

9{dU){p) < (respectively 9{dU){p) > 0) . 

This coincides with the usual inequality 6* < (0 > 0) on any subset of dD which 
is C^. Finally 9{dD) = in the viscosity sense if for any point p G dD we have 
both 9{dD) < in the viscosity sense, and 9{dD) > in the viscosity sense. 

It is tempting to conjecture the following: if fl is the set covered by smooth 
compact S"s satisfying 

9{S) := A + {g^^ - n'n^)K,j < , 
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then any outermost connected component of dQ is smooth, separating, and 
satisfies (1.22). When Kij = this statement is known to be true by the aheady 



mentioned result of G. Huisken and T. Ilmanen. 

We refer the reader to |^,|7^ and references therein for a numerical treatment 
of apparent horizons. 

1.7 Classification of stationary solutions ("No hair 
theorems") 

Uniqueness theorems for stationary solutions have been one of the central fields 
of research in the theory of black holes, the reader is referred to ^,34| and 



references therein for a review, as well as lists of open problems; cf. also [ |75| , 
Actually the perspective here is somewhat larger, as one is interested in classify- 
ing all solutions, not only those that contain black holes. Since the time of writing 
of ||3^,Q some progress has been made, both for zero and non-zero cosmological 
constant; it is convenient to discuss various cases separately. Throughout this 
section, stationary means that there exists an action of R on ^ by isometrics, 
with the orbits of the group being timelike sufficiently far away in an asymp- 
totic region; strictly stationary means that the orbits are timelike everywhere. 
Stationary black hole solutions are never strictly stationary. 

yl = 0: The following fundamental theorem has been recently proved by M. An- 
derson: 

Theorem 3 (Anderson ||^) The only geodesically complete strictly stationary 
solutions of the vacuum Einstein equations which do not contain closed timelike 
curves are the Minkowski space-time and quotients thereof. 

This is a milestone result in general relativity, involving only the natural 
geometric conditions of completeness together with absence of closed timelike 
curves. This should be contrasted with the version of this result originating in 
Lichnerowicz's work, the state-of-the art form of which, as based on previous 
techniques, assumes in addition global hyperbolicity, and asymptotic flatness, 
and existence of a compact interior. Last but not least, the paper introduces 
techniques which have not been used in mathematical general relativity so far, 
and which are likely to be useful tools in future work. 

In the black hole case, the only significant progress made since the review 
paper is for static black holes: In |^ the classification of static vacuum black 
holes which contain an asymptotically flat spacelike hypersurface with compact 
interior has been, essentially, finished. This problem has a long history, start- 
ing with the pioneering work of Israel [pl|] . The most complete result existing 

Some comments about the qualification "essentially" are in order: uniqueness theo- 
rems usually proceed in two steps, the first being the reduction of the problem to a 
Riemannian one, the second consisting of a uniqueness theorem for the Riemannian 
problem. The latter part of the proof in jsoj seems to be optimal; one could consider 
improving the former, see pa for a detailed discussion. 
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in the literature prior to |^6| was that of Bunting and Masood-ul-Alam 
who showed, roughly speaking, that all appropriately regular such black holes 
which do not contain degenerate horizons belong to the Schwarzschild family. 
In the condition of non-degeneracy of the event horizon is removed, show- 
ing that the Schwarzschild black holes exhaust the family of all appropriately 
regular black hole space-times. The proof requires an extension of the positive 
mass theorem which applies to asymptotically flat complete Riemannian mani- 
folds, proved in |ll|| . The paper contains an improvement of similar previous 
results concerning the electro-vacuum black holes, under the restrictive condi- 
tion that all degenerate components of the black hole carry charges of the same 
sign. This seems a reasonable condition from a physical point of view: oppo- 
site charges would attract, which added to the gravitational attraction should 
prevent the configuration from being static. However, a proper mathematical 
treatment which would transform this argument into a proof is still missing. 



A > 0: The natural boundary conditions for stationary solutions with yl > are 
the "no boundary ones" : one considers globally hyperbolic space-times contain- 
ing compact, boundaryless spatial surfaces. In this case there is no asymptotic 
region, and stationarity is defined by the requirement that the set of points at 
which the Killing vector field is timelike is non-empty. In | |9^ J. Lafontaine and 
L. Rozoy have announced a complete classification of static space-times with 
a compact (boundaryless) totally geodesic spacelike hypersurface S^, under the 
assumption that the Lorentzian norm ^/-gf,^X''X'' of the Killing vector field 
is a Morse-Bott function, and that of analyticity^ of the metric on S^. The 
hypothesis that =5^ is totally geodesic is the usual hypothesis of staticity. This 
is the first result of such generality in this context. The proof is an extension 
of an argument of R. Beig and W. Simon ||l^] done in a related context, and 
proceeds by proving that the metric on .y must be conformally flat; one can 
then conclude using ||9ll] . 

Uniqueness results under the completely different hypothesis of existence of 
a conformal completion have been previously established by H. Friedrich | |5^ 
and by Boucher . 



yl < 0: Our understanding of the class of stationary space-times with a negative 
cosmological constant is much poorer than that of the A = ones. It seems 
that the only results available so far concern static space-times. To start with, 
the question of boundary conditions which should be imposed in the asymptotic 
region are far from being understood, a beginning of a systematic analysis of this 
question can be found in |^5| . A uniqueness argument for the de Sitter metric 
in the strictly stationary case has been proposed in [H . It is only recently that 



complete proofs of the required positive energy theorem have been given |4lHl34 



This is a condition on the analyticity of the metric at the set of points at which 
the Killing vector vanishes or becomes null, since the metric is necessarily analytic 
elsewhere. 
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136 1; the results in pl| do actually lead to a uniqueness result under much weaker 
asymptotic conditions than originally suggested |jl^. It should be pointed out 
that the nature of the corresponding Riemannian problem is completely different 
from that of the asymptotically flat case: here one needs to analyze the set of 
solutions of the equations for a Riemannian metric g and a function y on a three 
dimensional manifold S'', 

AV = -AV , (1.24) 
% = V-'D,D,V + Ag,, , (1.25) 

and it is customary to assume that (J^, g) can be conformally compactified. Now, 
the uniqueness theorems mentioned require that V~^g extends by continuity to a 
smooth metric on the conformal boundary, with constant scalar curvature there. 
No reasons are known why this should be the case, and it is conceivable that 
more general static solutions would exist which do not satisfy this condition. It 
would be of interest to construct such solutions, or to prove their non-existence. 

The topology of the boundary of the conformally compactifiable solutions of 
(1.24)-( 1.25| ) does not have to be spherical, as is the case for the anti-de Sitter 



metric, while the uniqueness result mentioned above assumes that the boundary 
at infinity doo-^ is a sphere. When doo-^ is a torus, a strictly static solution 
has been found by Horowitz and Myers |Q . The metric in n + 1 > 4 spacetime 
dimensions is 

n-2 

= ^r'^dt^ + —-dr^ + V{r)dcb'' + Y , (1.26) 

V{r) ^ 

where Vir) = , = — '^^'Xa^ ■• ^"^^ p is a constant. Regularity de- 

mands that <\) be identified with period /3o — The periods of the ?/"s are 

arbitrary. The space-time metric is "asymptotically locally anti-de Sitter" with 
boundary at conformal infinity foliated by spacelike (n— l)-tori. The time slices of 
space-time itself, when conformally completed, have topology x T"~^, which 
in dimension 3-1-1 is a solid torus. Under a convexity condition on the confor- 
mal boundary, G. Galloway, S. Surya and E. Woolgar have proved a uniqueness 
theorem for the corresponding metrics [ ]62| . 

The above mentioned results concern the strictly static case, and the ques- 
tion of classification of black holes with a Killing vector which is not timelike 
everywhere is essentially open. In |45| ] an argument has been presented which 
shifts the problem of proving uniqueness of a class of such black holes to that 
of proving a Penrose-type inequalityP| for conformally compactifiable Rieman- 
nian manifolds with minimal boundary, and with scalar curvature bounded from 
below by a negative constant. While there is some hope that some such results 
could be proved by extending the arguments of [pO|, this remains to be seen. 



A version of such an inequality has been oriKinally proposed by Gibbons j^^; the 
inequality proposed there is not quite correct, see |45||. 
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1.8 Black holes without Scri 



There has been considerable progress in the numerical analysis of black hole 
solutions of Einstein's equations; here one of the objectives is to write a stable 
code which would solve the full four dimensional Einstein equations, with initial 
data containing a non-connected black hole region that eventually merges into a 
connected one. One wishes to be able to consider initial data which do not possess 
any symmetries, and which have various parameters - such as the masses of the 
individual black holes, their angular momenta, as well as distances between them 
- which can be varied in significant ranges. Finally one wishes the code to run to 
a stage where the solution settles to a state close to equilibrium. The challenge 
then is to calculate the gravitational wave forms for each set of parameters, 
which could then be used in the gravitational wave observatories to determine 
the parameters of the collapsing black holes out of the observations made. This 
program has been being undertaken for years by several groups of researchers, 
with steady progress being made |l|,|,||,|l|,[l9[|^,|7l|,|6|,||.0 

There is a fundamental difficulty above, of deciding whether or n ot on e is 
dealing indeed with the desired black hole initial data: the definition ( 1.13|) of 



a black hole requires a conformal boundary ^ satisfying some - if not all 



properties discussed in Sect. |l.4.1| . Clearly there is no way of ensuring those 
requirements in a calculation performed on a finite space-time grid.0 

In practice what one does is to set up initial data on a finite grid so that 
the region near the boundary is close to flat (in the conformal approach the 
whole asymptotically flat region is covered by the numerical grid, and does not 
need to be near the boundary of the numerical grid; this distinction does not 
affect the discussion here). Then one evolves the initial data as long as the 
code allows. The gravitational waves emitted by the system are then extracted 
out of the metric near the boundary of the grid. Now, our understanding of 
energy emitted by gravitational radiation is essentially based on an analysis of 
the metric in an asymptotic region where g is nearly flat. In order to recover 
useful information out of the numerical data it is thus necessary for the metric 
near the boundary of the grid to remain close to a flat one. If we want to be sure 
that the information extracted contains all the essential dynamical information 
about the system, the metric near the boundary of the grid should quiet down to 
an almost stationary state as time evolves. Now, it is straightforward to set-up 



Some spectacular visualizations of the calculations performed can be found at the 



URL jittp : // j ean-luc . aei . mpg . de/NCSA1999/GrazingBlackHoles 

The conformal approach developed by Friedrich (cf., e.g., jscj] and references therein, 
as well as S. Husa's and J. Frauendiener's contributions to this volume) provides an 
ideal numerical framework for studying gravitational radiation in situations where 
the extended space-time is smoothly conformally compactifiable across i"*", since then 
one can hope that the code will be able to "calculate Scri" globally to the future of 
the initial hyperboloidal hypersurface. It is not clear whether a conformal approach 
could provide more information than the non-conformal ones when is itself a 
singularity of the conformally rescaled equations, as is the case for black holes. 
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a mathematical framework to describe such situations without having to invoke 
conformal completions, this is done in the next section. 



1.8.1 Naive black holes 

Consider a globally hyperbolic space-time ^ which contains a region covered by 
coordinates {t, a;') with ranges 



r := ^^0, To-RQ + r<t<(X), (1.27) 

such that the metric g satisfies there 

< Cir-" < C2 , a>0, (1.28) 

for some positive constants Ci, C2, a; clearly C2 can be chosen to be less than 
or equal to CiRq". Making Rq larger one can thus make C2 as small as desired, 
e.g. 

C2 - 10-2 , (1.29) 

which is a convenient number in dimension 3 + 1 to guarantee that o bjects 
algebraically constructed out of g (such as g'^'^ , y/det g) are well controlled; ( 1.29|) 



is certainly not optimal, and any other number suitable for the purposes at hand 
would do. To be able to prove theorems about such space-times one would need 
to impose some further, perhaps not necessarily uniform, decay conditions on a 
finite number of derivatives of g; there are various possibilities here, but we shall 
ignore this for the moment. Then one can define the exterior region ^ext, the 
black hole region and the future event horizon S" as 



^ext 



:=Ur>ToJ {Sr,R„) = J {Dr>T„Sr,Ro) , (1-30) 

^:=^\^oxt, <^:=d£§, (1.31) 

where 

Sr.Ro-={t = r,r^RQ} . (1.32) 



We will refer to the definition ( 1.27 )-( 1.32 ) as that of a naive black hole 



In the setup of ( 1.27 )-( p".32| ) an arbitrarily chosen Rq has been used; for this 
definition to make sense ^ so defined should not depend upon this choice. This 
is indeed the case, as can be seen as follows: 

Proposition 4 Let C \ 5(0, i?o), a = 1,2, and let X C ^ be of the 
form {(t > Tq — RQ + r{x),x) ,x G i^q} in the coordinate system of (1.28). Then 
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Proof. If 7^ is a future directed causal path from p E ^ to g = (i, a;) e "^i, then 
the path obtained by concatenating F with the path [0, 1] 9 s — > (t(s) := i + 
s, a;(s) := x) is a causal path which is not a null geodesic, hence can be deformed 
to a timelike path from p to + 1, cc) e "^i. It follows that 7~('^i) = 
clearly the same holds for '^2. Next, let G ^a, and let 7 : [0, 1] ^ R3\B(0, i?o) 
be any differentiable path such that 7(0) = X\ and 7(1) — x^. Then for any > 
To — Ro + r{xi) the causal curve [ 0, 1] 3 s r{s) = {t := Cs + to,x{s) :— 7(5)) 
will be causal for the metric g by ( 1.28| ) if the constant C is chosen large enough, 
with a similar result holding when xi is interchanged with X2- The equality 
J~('^i) = I~{'^2) easily follows from this observation. □ 



Summarizing, Prop. ^ shows that there are many possible equivalent definitions 
of ^e^t- in ( 1.30| ) one can replace J^{Sr,R„) by J~{Sr,Ri) for any Ri > Rq, but 
also simply by J^{{t + r, q)), for any p — {t, q) G ^ which belongs to the region 
covered by the coordinate system {t,x''). 

The following remarks concerning the definition (1.3C)-(1.31) are in order: 



For vacuum, stationary, asymptotically flat space-times the definition is 
equivalent to the usual one with Footnote 7, p. 572]; here the re- 

sults of [^|5^ are used. H owev er , one does not expect the existence of a 
smooth ^+ to follow from ( 1.27 )-( 1.28 ) in general. 

Suppose that ^ adm its a conformal completion in the sense defined at the 
beginning of Sect. lA , and that .y is semi- complete to the fut ure, in the sense 
that the Geroch-Horowitz condition of the beginning of Sect. 1.4.1 holds with 
y ^ M. X there replaced by D M+ x S^. Then for any finit e int erval 
[To, Ti] there exists Ro{Tq, Ti) and a coordinate system satisfying ( 1.28| ) and 
covering a set r > Ro{To,Ti), Tq - Rq < t < Ti - Rq. This follows from 
the Tamburino-Winicour construction of Bondi coordinates (w, r, 0, ip) near 
^+ [124 , followed by the introduction of the usual Minkowskian coordinates 
t — u + r, X = r sinO cosip, etc. The problem is that R{Ti, T2) could shrink to 
zero as T2 goes to infinity. Thus, when ^+ exists, conditions (1.27)-( 1.28| ) are 
uniformity conditions on ^+ to the future - the metric remains uniformly 
controlled on a uniform neighbourhood of as the retarded time goes to 
infinity. 

It should not be too difficult to check whether or not the future geodesically 

and of Christodoulou and Klain- 
~2i, admit coordinate 



complete space-times of Friedrich [ |53| , 
erman [ |3l[ |, or the Robinson- Trautman black- holes 
systems satisfying ( 1.27D -(1.28). 



It is not clear if asymptotically flat space-times in which no such control is 
available do exist at all; in fact, it is tempting to formulate the following version 
of the (weak) cosmic censorship conjecture: 
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The maximal globally hyperbolic development of generic^, asymptoti- 
cally flat, vacuum initial data contains a region with coordinates satis- 



tying (1.27)-(1.28) 



Whatever the status of this conjecture, one can hardly envisage numerical sim- 
ulations leading to the calculation of an essential fraction of the total energy 
radiated away in space-times in which some uniformity conditions do not hold. 



1.8.2 Quasi-local black holes 

As already argued, the naive approach of the previous section should be more 
convenient for numerical simulations of black hole space-times, as compared to 
the usual one based on Scri. It appears to be even more convenient to have a 
framework in which all the issues are localized in space; we wish to suggest such 
a framework here. When numerically modeling an asymptotically flat space- 
time, whether in a conformal or a direct approach, a typical numerical grid 
will contain large spheres S(R) on which the metric is nearly flat, so that an 
inequality such as (1.28) -([l.29[ ) will hold in a neighbourhood of S(R). On slices 



t — const the condition ( |l.28| ) is usually complemented with a fall-off condition 
on the derivatives of the metric 

< C'^"""' . (1-33) 



however condition ( 1.33 ) is inadequate in the radiation regime, where retarded 
time derivatives of the metric are not expected to fall-off faster than r~^. It 
turns out that there is a condition on derivatives of the metric in null directions 
which is fulfilled at large distance both in spacelike and in null directions: let 
Ka, a = 1,2, be null future pointing vector fields on S{R) orthogonal to S{R), 
with Ki inwards pointing and K2 - outwards pointing; these vector fields are 
unique up to scaling. Let da denote the associated null second fundamental forms 
defined as 

yX,YeTS{R) Xa{X,Y):^g{VxKa,Y) . (1.34) 

It can be checked, e.g. using the asymptotic expansions for the connection coef- 
ficients near ^+ from Appendix C], that xi is negative definite and X2 is 
positive definite for Bondi spheres S{R) sufficiently close to J^~^; similarly for 
J^^. This property is not affected by the rescaling freedom at hand. Following 
G. Galloway |^0|, a two-dimensional spacelike submanifold of a four-dimensional 
space-time will be called weakly null convex if xi is semi-positive definite, with 
the trace of X2 negative.|^ The null convexity condition is easily verified for 



sufficiently large spheres in a region asymptotically fiat in the sense of (1.33) 



The examples constructed by Christodoulou [^oj with spherically symmetric gravi- 
tating scalar fields suggest that the genericity condition is unavoidable, though no 
corresponding vacuum examples are known. 

Galloway defines null convexity through the requirement of positive definiteness of 
Xi and negative definiteness of X2 - However, he points out himself p. 1472] that 
the weak null convexity as defined above suffices for his arguments to go through. 
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It does also hold for large spheres in a large class of space-times with negative 
cosmological constant. The null convexity condition is then the condition which 
we propose as a starting point to defining "quasi-local" black holes and horizons. 
The point is that several of the usual properties of black holes carry over to the 
weakly null convex setting. In retrospect, the situation can be summarized as 
follows: the usual theory of Scri based black holes exploits the existence of con- 
jugate points on appropriate null geodesies whenever those are complete to the 
future; this completeness is guaranteed by the fact that the conformal factor goes 
to zero at the conformal boundary at an appropriate rate. Galloway's discovery 
in is that weak null convexity of large spheres near Scri provides a second, 
in principle completely independent, mechanism to produce the needed focusing 
behaviour. 

Throughout this section we will consider a globally hyperbolic space-time 
(^, g) with time function t. Let =f7 C ./# be a finite union of connected timelike 
hypersurfaces 5^ in we set 

{t = r} , ^(r) ^ n , 5L(t) := X^HS^r ■ (1-35) 

For further purposes anything that happens on the exterior side of X is com- 
pletely irrelevant, so it is convenient to think of =^ as a boundary of global 
hyperbolicity should then be understood in the sense that :— ^ U ^,51) is 
strongly causal, and that J^{p; ^)r]J~{q; ^) is compact in ^ for allp, q S ^ . 
Recall that the null convergence condition is the requirement that 

Ric(X, X) > for aU X e TM . (1.36) 

We have the following topological censorship theorem for weakly null convex 
timelike boundaries: 



Theorem 5 (Galloway [|60[ ) Suppose that a globally hyperbolic space-time , g) 
satisfying the null convergence condition ( l.Sdj) has a timelike boundary 3' — 



^i,=i^a o,nd a time function t such that the level sets of t are Cauchy surfaces, 
with each section X{t) of 3^ being null convex. Then distinct 3'a 's cannot com- 
municate with each other: 

a^P j+(^„)n J-(^^) =0 . 

As is well known, topological censorship implies constraints on the topology: 

Theorem 6 (Galloway [|60[ ) Under the hypotheses of Theorem^ suppose fur- 
ther that the cross-sections =5^(t) of 3'a have spherical topology^ Then the 
a-domain of outer communication 

((^„)> := J+(^„)n J-(^„) (1.37) 

is simply connected. 

The reader is referred to and references therein for results without the hypoth- 
esis of spherical topology. The results there, presented in a Scri setting, generalize 
immediately to the weakly null convex one. 
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It follows in particular from Theorem ^ that ^ can be replaced by a subset 
thereof such that ^ is connected in the new space-time, with all essential prop- 
erties relevant for the discussion in the remainder of this section being unaffected 
by that replacement. We shall not do that, to avoid a lengthy discussion of which 
properties are relevant and which are not, but the reader should keep in mind 
that the hypothesis of connectedness of 3" can indeed be done without any loss 
of generality for most purposes. 

We define the quasi-local black hole region and the quasi-local event 

horizon Sa;-^ associated with the hypersurface ^ by 



'\J-{X,) , Sg-^ :=9^^„ . (1.38) 



// 9' is the hypersurface U^>t'o5t-^;jo of Section 1.8.1 then the resulting black 



hole region coincides with that defined in (1.31), hence does not depend upon 
the choice of Rq by Proposition^; however, S^sr^ might depend upon the chosen 
family of observers ^ in general. It is certainly necessary to impose some further 
conditions on to reduce this dependence. A possible condition, suggested by 
the geometry of the large coordinate spheres considered in the previous section, 
could be that the light-cones of the induced metric on <f7 are uniformly controlled 
both from outside and inside by those of two static, future complete reference 
metrics on . However, neither the results above, nor the results that follow, 
do require that condition. 

The Scri-equivalents of Theorem | ||2|, |6[ ||, |6l|, ||, ||l allow one to control 
the topology of "good" sections of the horizon, and for the standard stationary 
black-holes this does lead to the usual S*^ x M topology of the horizon [^,^. In 
particular, in stationary, asymptotically flat, appropriately regular space-times 
the intersection of a partial Cauchy hypersurface with an event horizon will 
necessary be a finite union of spheres. In general space-times such intersections 
do not even need to be manifolds: for example, in the usual spherically symmetric 
collapsing star the intersection of the event horizon with level sets of a time 
function will be a point at the time of appearance of the event horizon. We refer 
the reader to Section 3] for other such examples, including one in which the 
topology of sections of horizon changes type from toroidal to spherical as time 
evolves. This behaviour can be traced back to the existence of past end points 
of the generators of the horizon. Nevertheless, some sections of the horizon have 
controlled topology - for instance, we have the following: 

Theorem 7 Under the hypotheses of Theorem consider a connected compo- 
nent 3'a of £^ such that Sg-^ ^ . Let 

^^{t) :=aj+(^„(r)) . 

if '^a (j) n Sy^ is a topological manifold, then each connected component thereof 
has spherical topology. 

Proof. Consider the open subset of ^ defined as 

:=/+C^„(T);^)n/-(^„;^) C . 
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We claim that (^^riSl^^) is globally hyperbohc: indeed, let p,q £ ^t', global 
hyperboHcity of ^ shows that J~ {p; Ji') n J+ (g; ^) is a compact subset of 
which is easily seen to be included in J^t. It follows that ^T)^J^{q\ J^t^ 

is compact, as desired. By the usual decomposition we thus have 

w R X 5 , 

where 5 is a Cauchy hypersurface for Applying Theorem ^ to the globally 
hyperbolic space-time (which has a weakly null convex boundary n {t > 
r}) one finds that is simply connected, and thus so is S. Since '^air) and 
Sa are null hypersurfaces in it is easily seen that the closure in ^ of the 
Cauchy surface {0} x S intersects precisely at '^a(r) n Sg-^. It follows that 
iS is a compact, simply connected, three dimensional topological manifold with 
boundary, and a classical result j74[ Lemma 4.9] shows that each connected 
component of dS is a sphere. The result follows now from dS ~ '^a (t) n Ssr^ . 
□ 

Yet another class of "good sections" oi Sg can be characterized^ as follows: 
suppose that ((=5^)) n S^r is a submanifold with boundary of ^ which is, more- 
over, a retract of ((.5^)). Then ((=3^)) n is simply connected by Theorem ^ 
and spherical topology of all boundary components of ((5^)) n S^r follows again 
from ItJ, Lemma 4.9]. It is not clear whether there always exist time functions 
t such that the retract condition is satisfied; similarly it is not clear that there 
always exist r's for which the conditions of Theorem ^ are met for metrics which 
are not stationary (one would actually want "a lot of r's"). It would be of interest 
to understand this better. 

We have an area theorem for S 3? : 

Theorem 8 Under the hypotheses of Theorem suppose further that Sg- ^ 
0. Let 5^ a, a = 1,2 he two achronal spacelike embedded hypersurfaces of 
differentiability class, set Sa — H (ag. Then: 

1. The area of Sa is well defined. 

2. If 

Si c J- (52) , 

then the area of S2 is larger than or equal to that of Si. (Moreover, this 
is true even if the area of Si is counted with multipliciti^ of generators 
provided that SiD S2 — 

We note that point 1 is less trivial as it appears, because horizons can be 
rather rough sets, and it requires a certain amount of work to establish that 
claim. 

I am grateful to G. Galloway for useful discussions concerning this question, as well 
as many other points presented in this section. 
'^^ See ]|| for details. 
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Proof. The result is obtained by a mixture of methods of and of |6^, and 
proceeds by contradiction: assume that the Alexandrov divergence 9m oi S"^ is 
negative, and consider the Se,ri,5 deformation of the horizon as constructed in 
Proposition 4.1 of with parameters chosen so that 6'^ a < 0. Global hyper- 
bolicity implies the existence of an achronal null geodesic from Se,ri,s to some cut 
■^5'(r) of c^'. The geodesic can further be chosen to be "extremal", in the sense 
that it meets ^{t) for the smallest possible value of t among all generators of 
the boundary of J~^{Se^n,5) meeting The argument of the proof of Theorem 1 
of |]60| shows that this is incompatible with the null energy condition and with 
weak null convexity of =^(t). It follows that 9^ > 0, and the result follows from 
Proposition 3.3 and Theorem 6.1]. □ 

It immediately follows from the proof above that, under the hypotheses of The- 
orem § the occurrence of twice differentiable future trapped (compact) surfaces 
implies the presence of a black hole region. The same result holds for semi-convex 



compact surfaces which are trapped in an Alexandrov sense, that is, (1.21) holds 
with A there defined in a way which should be clear from the discussion follow- 
ing Theorem |^. It is, however, not known if the existence of marginally trapped 
surfaces - whether defined in a classical, or Alexandrov, or a viscosity sense - 
does signal the occurrence of black hole; it would be of interest to settle that. 

In summary, we have shown that the quasi-local black holes, defined using 
weakly null convex timelike hypersurfaces, or boundaries, possess several prop- 
erties usually associated with the Scri-based black holes, without the associated 
problems. We believe they provide a reasonable alternative, well suited for nu- 
merical calculations. In the next section we address the question of numerically 
finding the resulting horizons. 



1.8.3 Finding horizons 



Consider the large coordinate spheres Sr,Ri of Section 1.8.1; under reasonable 



restrictions on ^ one would expect that dJ (iST,iJi) converges as r ^ oo — e.g. 



in Hausdorff topology — to the event horizon of (1.31). Similarly, the achronal 
boundaries dJ~ {£^{t)) should sometimes converge to the event horizon S^. 
Whenever that is the case, the 9J^(iS,-._Ri)'s or the (.5^(r))'s can be taken 
as good approximations of the event horizon when r is chosen to be large enough. 
This is of practical significance, as it provides a way of approximately locating 
the horizon in numerical simulations. Let us establish one such approximation 
result, for the (9J~(,5^(t))'s; clearly the corresponding result for the J~ [Sr.RiYs 
follows by specialization. We assume that ^ is spatially compact and globally 
hyperbolic, which seem rather natural assumptions in the context of numerical 
simulations of space-times: 

Theorem 9 Consider a globally hyperbolic space-time (.-#, g) with compact Cauchy 
hypersurfaces S^r and timelike boundary 3' , and suppose that 



(1.39) 
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Then for any t we have 

where the convergence is in HausdorjJ distance. 

Proof. Under the usual identification of ^ with R x .y^, let / be the graphing 
function of dJ~ over the projection pr2t^ of on S^q; prjtc? is an open 

subset of c5^o by the invariance-of-the-domain theorem. If p = {t,x) is such that 
t < f{x), then there exists a causal curve from p to 3^, hence p E J^{,^{a)) for 
some a. It follows that 

pr2(f -U^pr2J-(5^(a)) . 

Let /cr be the graphing function of dJ^{^{a)) over pr2 J^(=5^(cr)) C S'q; clearly 

pr2J-(^(a)) Cpr2J-(^(a')) a < a' , 

which shows that / and all the /o-/ 's are defined on pr2 J~ (>^(cr)) for a' > a. The 
/r's are uniformly Lipschitz and bounded over any compact subset K C V"^^^ of 
=5^0: E^nd hence for any sequence there exists a subsequence such that the /t-^ 's 
converge to a Lipschitz function h. Since fr 'Si f for all r we have h < f\K- Sup- 
pose that there exists x such that h{x) < f{x), then p = {{h{x) + f{x)) /2, x) 
is in J~{^), hence in J^{3^{a)) for some a, so that {h{x) + f{x)) /2 < fa{x) 
which is not possible; we thus have h = f\K- It follows that the /r's converge 
pointwise to /. The /^'s are monotonously increasing as r increases, and Dini's 
theorem implies that the /t-'s converge uniformly to / over any compact set. The 
result follows now by elementary considerations using the fact that S"^^ C\ S^r is 
compact, being the intersection of a closed set with a compact one. □ 

Theorem | gives partial justificat ion for the numerical analysis of Aninos 0, 



Libson et al. ||9^, or Masso et al. |102|,^ where the event horizon § is taken to 
be 

S K, dJ~ {St,Ri ) , T as big as the computer simulation allows . 

In those works dJ~{Sr,Ri) is further numerically approximated as the solution 
of an eikonal equation, which leads to difficulties at past end points of the gen- 
erators of dJ^ {St,Ri). We conjecture that a more stable method of locating 
objects such as dJ^ {St,f/.i) is provided by the following straightforward modifi- 
cati on of the Fermat-type variational principle presented in the Appen dix below, 
see drU): let V{x) be defined as in the paragraph preceding ( |l.43[) , with 



there replaced by Sr,Ri , and let 7 be the null lifts of paths 7 as defined in the 

We note that the claim in |^^, that the horizon acts as an attractor for backwards- 
directed null geodesies, does not seem to be justified. In fact, such a statement is 
coordinate-dependent: it is easy to introduce coordinate systems in the Schwarzschild 
space-time where backwards-directed null geodesies will be actually repelled by the 
horizon. 
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paragraph preceding ( 1.45| ). An argument similar to that of the proof of Corol- 



lary A. 3 shows that in globally hyperbolic space-times the set dJ {Sr.Ri) is a 



graph of a function / such that 

f(x) = sup i(7(a)) . (1.40) 

Such a maximization over null lifts 7 of paths 7 which have image in a fixed 
spacelike hypersurface should be easy to perform numerically by e.g. Newton's 
method. The variational principle above automatically takes care of the past 
end-points of the horizon. 

Acknowledgements The author is grateful to Y. Choquet-Bruhat, G. Horowitz, 
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improvements to a previous version of this manuscript, and to R. Geroch and 
J. Jezierski for useful comments and discussions. 



Appendix A: Future horizons are semi-convex 

A hypersurface C M will be said to be future null geodesically ruled if ev- 
ery point p £ Jif belongs to a future inextensible null geodesic F C J^; those 
geodesies will be called the generators of Jf. We emphasize that the genera- 
tors are allowed to have past endpoints on Jif, but no future endpoints. Past 
null geodesically ruled hypersurfaces are defined by changing the time orienta- 
tion. Examples of future geodesically ruled hypersurfaces include past Cauchy 



horizons V (E) of achronal sets S [115, Theorem 5.12] and black hole event 
horizons dJ-{y+) @ p. 312]. 

Let dimM = n + 1 and suppose that O is a domain in R". Recall that 
a continuous function / : O — > R is called semi-convex if there exists a 
function : O ^ R such that / + is convex. We shall say that the graph 
of / is a semi-convex hypersurface if / is semi-convex. A hypersurface in a 
manifold M will be said semi-convex if can be covered by coordinate patches 
Ua such that Jif (iLia is a semi-convex graph for each a. 

Consider an achronal hypersurface 7^ in a globally hyperbolic space- 
time (M, g) . Let t be a time function on M which induces a diffeomorphism of M 
with Rx E in the standard way [|4, 122|, with the level sets Et = {p\t{p) — r} of 



t being Cauchy surfaces. As usual we identify Eq with E, and in the identification 
above the curves R x {q}, q £ E, are integral curves of Vt. Define 

E^ = {g e ^|R X {q} intersects Jf"} . (1.41) 

For q E Ejif the set {I x q) O H is a point by achronality of Jf, we shall denote 
this point by (/(<?), q). Thus an achronal hypersurface in a globally hyperbolic 
space-time is a graph over of a function /. The invariance-of-the-domain 
theorem shows that Ejfc is an open subset of E. We have the following: 
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Theorem A.l Let .jZ/f ^ ^ be an achronal future null geodesically ruled hyper- 
surface in a globally hyperbolic space-time (M = Mx 17, (7). Then is the graph 
of a semi-convex function f defined on an open subset of S, in particular 
is semi-convex. 

Proof. As discussed above, is the graph of a function /. The idea of the proof 
is to show that / satisfies a variational principle, the semi-concavity of / follows 
then by a standard argument. Let p g .J^ and let C be a coordinate patch in a 
neighborhood of p such that x'^ — t, with O of the form / x B{3R), where B{R) 
denotes a coordinate ball centered at of radius R in R^, with p ~ {t{p),0). 
Here / is the range of the coordinate a;*^, we require it to be a bounded interval 
the size of which will be determined later on. We further assume that the curves 
/ X {x}, X E B{3R), are integral curves of Vi. Define 

Ua = {x E B{3R)\ the causal path I 3 t ^ {t,x) intersects J^} . 

We note that Uq is non-empty, since EUq. Set 

J^^=J^nSa. (1.42) 

and choose a large enough so that O C I~'{X!a)- Now p lies on a future inexten- 
sible generator F of J^f, and global hyperbolicity of {M, g) implies that F n 2J„ 
is nonempty, hence Jifa is nonempty. 

For x E B{3R) let V{x) denote the collection of piecewise differentiable 
future directed null curves F : [a, b] M with F{a) e M x {x} and F{b) E M'a- 
We define 

t{x) = sup t{F(a)) . (1.43) 
We emphasize that we allow the domain of definition [a, 5] to depend upon F , 



and that the "a" occurring in t(F{a)) in (1.43) is the lower bound for the domain 



of definition of the curve F under consideration. We have the following result 
(compare ||, |8 , 116 1 ) : 



Proposition A. 2 (Fermat principle) For x E Uq we have 

t{x) = f{x) . 

Proof. Let F be any generator of J^f such that F{0) = {f{x),x), clearly F E 
V{x) so that t{x) > f{x). To show that this inequality has to be an equality, 
suppose for contradiction that t{x) > f{x), thus there exists a null future di- 
rected curve F such that t{F{0)) > f{x) and F{1) E Jf^ C . Then the curve 
F obtained by following M x {a;} from {f{x),x) to {t{F{0)),x) and following F 
from there on is a causal curve with endpoints on Jif which is not a null geodesic. 
By Proposition 4.5.10] F can be deformed to a timelike curve with the same 
endpoints, which is impossible by achronality of J^. □ 



The Fermat princ iple. Proposition A. 2, shows that / is a solution of the vari- 
ational principle (1.43). Now this variational principle can be rewritten in a 
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somewhat more convenient form as follows: The identification of M with M x 17 
by flowing from = S along the gradient of t leads to a global decomposition 
of the metric of the form 

g = a{-df + ht) , 

where ht denotes a t-dependent family of Riemannian metrics on S. Any future 
directed differentiable null curve r{s) = {t{s),^{s)) satisfies 



^ = v/'^t(.)(7,7), 



ds 

where 7 is a shorthand for d^{s)/ds. It follows that for any F £ 'P{x) it holds 
that 



ma)) = tir{b))- f^ds 

Ja ds 

\Jht(s){i,i)ds . 



This allows us to rewrite (1.43) as 

t[x) = cr - ii{x) , fi{x) = ^M^^ J ^ht(s){i,i)ds . (1.44) 

We note that in static space^imes ii{x) is the Riemannian distance from x to 
J^a-- In particular Equation (1.44) implies the well known fact, that in globally 
hyperbolic static space-times Cauchy horizons of open subsets of level sets of t 
are graphs of the distance function from the boundary of those sets. 

Let 7 : [a, 6] — > be a piecewise differentiable path, for any p e Mx {7(6)} we 
can find a null future directed curve ^ : [a,b] M of the form 7(5) = (0(s), 7(3)) 
with future end point p by solving the problem 

r </)(&) = tip) , 

The path 7 will be called the null lift of 7 with endpoint p. 

As an example of application of Proposition A. 2 we recover the following well 
known result [115|: 

Corollary A. 3 / is Lipschitz continuous on any compact subset of its domain 
of definition. 

Proof. For y,z E B{2R) let C M x B{2R) be a compact set which contains all 
the null lifts Fy^z of the coordinate segments [y, z] := {\y+{l — X)z , A e [0, 1]} 
with endpoints {t{z)^z). Define 



C ^sn^{Jhp{n,n)\peK,\n\s = 1} , (1.46) 
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where the suprcmum is taken over all points p Cz K and over all vectors n G TpM 
the coordinate components of which have Euclidean length \n\s equal to one. 
Choose / to be a bounded interval large enough so that K C I x B{2R) and, as 
before, choose a large enough so that I x B{2R) lies to the past of Ua- Let y,z e 
B{2R) and consider the causal curve F = (i(s),7(s)) obtained by following the 
null lift Fy z in the parameter interval s G [0, 1], and then a generator of 
from {t{z),z) until in the parameter interval s E [1,2]. Then we have 



Further F eVix) so that 



m(2/) < J \lht{s){i,i)ds 

= \Jht{s){i.i)ds + \Jht(s){i,i)ds 



< C\y - z\s + ^JL{z) , (1.47) 



where | • \s denotes the Euclidean norm ofa vector, and with C defined in ( 1.46| ). 



Setting 1) first y = x, z = x + h in ( 1.47 ) and 2) then z = x, y = x -\- h, the 



Lipschitz continuity of / on B{2R) follows. The general result is obtained now 
by a standard covering argument. □ 

Returning to the proof of Theorem [A.l| , for x e B{R) let be a generator of 
such that Fx{Q) — {t{x),x), and, if we write Fx{s) — (4>x{s),jx{s)), then 
we require that -fx{s) e B{2R) for s e [0,1]. For s e [0,1] and h e B{R) let 
lx,±{s) € E be defined by 

7x,±(s) - lx{s) ± (1 - s)h = s^x{s) + (1 - s)(7,(s) ± /i) e B{2R) . 



We note that 



7x,±(0) = a; ± , 7x,±(l) = 7x(l) , 7x,± - 7x = T/j- • 

Let -Ck,± = {,4'±{s),lx,±) be the null lifts of the paths 7a;, ± with endpoints 
rx{l). Let if be a compact set containing all the Ilc,±'s, where x and h run 
through B{R). Let I be any bounded interval such that / x B{2R) contains K. 
As before, choose a so that / x B{2R) lies to the past of Sa, and let b be such 
that Fx{h) € J^a- (The value of the parameter h will of course depend upon 
x). Let F± be the null curve obtained by following Fx^± for parameter values 
s £ [0, 1], and then Fx for parameter values s G [1, h]. Then F± E V{x ± h) so 
that we have 

li{x ±h)< J y^ft.0±(s)(7n;,±7 7x,±)rfs + ^ \Jh^^(s)['^x,ix)ds . 



36 P. T. Chrusciel 



Further 

hence 

fi{x + h) + fi{x — h) 



^ \/h.4,^(s)iix,ix) I ds . 

(1.48) 



Since solutions of ODE's with parameters are differentiable functions of those, 
we can write 



</.±(s) = <^,(s) + 4,,{s)h' + ris, h), |r(s, h)\ < C\h\1 



(1.49) 



for some functions ipi, with a constant C which is independent of x,h d B{R) 
and s E [0,1]. Inserting ( L.49| ) in ( 1.48 ), second order Taylor expanding the 

function ^Jh^^(^g^{jx,±,ix,±){s) in all its arguments around (0a,(s), 7a,(s), 7a,(s)) 

and using compactness of K one obtains 



^{x + h) + fj,{x — h) 



tx{x)<C\h\l , 



(1.50) 



for some constant C. Set 



^{x) = ^i{x) - C\x\l 



Equation (1.50) shows that 



\/x,heB{R) ip{x)> 



ip{x + h) + tpix - h) 



A standard argument implies that ip is concave. It follows that 

f{x) + C\x\l - t{x) + C\x\l - a - /i(a;) + C\x\l = a - 4>{x) 
is convex, which is what had to be established. □ 
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